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ABSTRACT 


A computational investigation of the aerodynamic effects on fluidic thrust vectoring has been conducted. Three- 
dimensional simulations of a two-dimensional, convergent-divergent (2DCD) nozzle with fluidic injection for pitch 
vector control were run with the computational fluid dynamics code PAB using turbulence closure and linear 
Reynolds stress modeling. Simulations were computed with static freestream conditions (A/=0.05) and at Mach 
numbers from A/=0. 3 to 1.2, with scheduled nozzle pressure ratios (from 3.6 to 7.2) and secondary to primary total 
pressure ratios of p ( /p tp =0.6 and 1 .0. Results indicate that the freestream flow decreases vectoring performance and 
thrust efficiency compared with static (wind-off) conditions. The aerodynamic penalty to thrust vector angle ranged 
from 1.5 degrees at a nozzle pressure ratio of 6 with Af=0.9 freestream conditions to 2.9 degrees at a nozzle pressure 
ratio of 5.2 with M=0.7 freestream conditions, compared to the same nozzle pressure ratios with static freestream 
conditions. The aerodynamic penalty to thrust ratio decreased from 4 percent to 0.8 percent as nozzle pressure ratio 
increased from 3.6 to 7.2. As expected, the freestream flow had little influence on discharge coefficient. 


INTRODUCTION 

There are a host of potential benefits for 
incorporating thrust vectoring into military aircraft: 
enhanced turn rate, improved maneuverability, vertical 
and short take-off and landing (V/STOL) capabilities, 
and elimination of control surfaces, to name a few. 
Eliminating conventional control surfaces in lieu of 
propulsive controls may result in weight, signature, 
and drag reductions. The use of fluidic injection for 
thrust vectoring, instead of mechanized hardware, offers 
further improvements in aircraft signature, weight, and 
drag. 

The thrust vectoring studies initiated in the 1960's 
included gimbal and swivel mechanisms added to the 
tailpipe, prior to the axisymmetric nozzle, to improve 
take-off and landing performance (refs. 1-3). These 
mechanisms proved efficient for thrust vectoring under 
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static conditions, but resulted in detrimental penalties 
for installed conditions. 

The studies throughout the 1970's investigated a 
variety of non-axisymmetric nozzles: such as two- 
dimensional convergent-divergent (2DCD), wedge, and 
single-expansion-ramp nozzles (SERN), for their 
ability to accommodate thrust vectoring (refs. 1-10). 
Non-axisymmetric, thrust-vectoring nozzles offered 
integration and installed performance benefits over the 
conventional axisymmetric nozzles (ref. 1 1). 

Mechanical, multiaxis thrust vectoring was 
investigated in the 1980’s. Augmenting conventional 
controls with pitch and yaw propulsive control allowed 
for supermaneuverability or operation in the post-stall 
regime, where conventional controls are ineffective 
(refs. 12-13). Pitch and yaw vectoring was 
accomplished primarily with deflection of the divergent 
flaps, but integration concepts and post-exit vanes 
were also investigated (refs. 14-18). Most concepts 
that vectored the flow with subsonic turning achieved 
the best vectoring effectiveness and incurred the least 
thrust performance penalties. Installed performance of 
mechanical multiaxis thrust vectoring nozzles was 
shown to be mostly independent of angle-of-attack and 
sideslip angle (ref. 19), unlike conventional 
aerodynamic control surfaces. 
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NOMENCLATURE 


Fluidic injection for thrust vectoring was 
investigated in the 1990's (ref. 20-26). Instead of 
deflections of divergent flaps or vanes to create 
pressure differentials and vectored thrust, fluidic thrust 
vectoring introduces a secondary air stream into the 
primary jet flow to create an off-axis deflection of the 
jet thrust. Two of the most mature fluidic techniques 
are throat skewing and shock vector control. However, 
some less mature concepts like counterflow (ref. 25- 
26) and synthetic jets (ref. 27) have also been 
introduced for thrust vectoring. Throat skewing was 
built on the database of mechanical concepts that 
demonstrated the efficiency of turning the flow 
subsonically. Fluidic injection occurs near the throat 
of the nozzle to shift the sonic line, turn the flow 
subsonically, and create an asymmetric pressure 
loading in the throat-skewing concept. For the shock 
vector control technique, the secondary air stream is 
introduced in the divergent portion of the nozzle and 
hence, affects supersonic flow. Similar to the 
mechanical concepts that turn the flow supersonically, 
the shock vector control technique has more thrust 
losses associated with it, but in general, substantial 
vector angles can be achieved. In addition, fluidic 
thrust vectoring is generally more effective at 
overexpanded nozzle pressure ratios. To date, most of 
the data on fluidic thrust vectoring has been obtained at 
static freestream conditions (refs. 20-26). Unlike 
mechanical thrust vectoring, there is very little, if any, 
data documenting the effects of the freestream flow on 
thrust vectoring with fluidic injection. 

The current investigation attempted to initiate a 
database of external flow effects on fluidic thrust 
vectoring. The nozzle under investigation was a two- 
dimensional, convergent-divergent (2DCD) nozzle 
with fluidic injection for pitch thrust vector control. 
The secondary air stream was injected through a slot 
in the lower divergent wall. Computational fluid 
dynamics (CFD) was used to simulate the nozzle at 
typical nozzle pressure ratio (NPR) schedule based on 
Mach number, up to M-\2. Simulations were also 
computed at the same NPR schedule with static 
freestream conditions (M=0.05) to document the effect 
of the external freestream on vectoring effectiveness, 
thrust efficiency, and discharge coefficient. Secondary 
to primary total pressure ratios of 0.6 and 1 .0 were 
investigated, which corresponded to secondary weight 
flow rates of 5.8 percent and 9.65 percent of the 
primary weight flow rate, respectively. 
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exit area, 13.104 in 2 
slot minimum area, 0.565 in 2 
throat area, 5.416 in 2 
expansion ratio, 2.4 

discharge coefficient, (w s +w p ) / ( w, v + w ip ) 
primary flow discharge coefficient, w p / w ip 
secondary flow discharge coefficient, w s / w, 
thrust ratio for no injection cases, F r / F ip 
system thrust ratio, F r / (F ip + F i s ) 
turbulence viscosity coefficient, 0.09 


axial force, lb 
normal force, lb 
resultant force, lb 

ideal isentropic thrust of primary flow, lb 

ideal isentropic thrust of secondary flow, lb 

turbulent kinetic energy 

aerodynamic body length, 64.745 inches 

free stream Mach number 

model station, inches 

unit normal vector, (n„ n 2 , n^) 

nozzle pressure ratio, p^j/p^ 

design nozzle pressure ratio, 14.588 

surface static pressure, psi 

average jet total pressure, psi 

primary flow total pressure, psi 

secondary flow total pressure, psi 

secondary to primary total pressure ratio 

free-stream static pressure, psi 

free-stream dynamic pressure, psi 

strain component, 1/sec 

velocity vector 

actual primary weight flow 

ideal primary weight flow 

actual secondary weight flow 

ideal secondary weight flow 

normalized axial location along afterbody 


pitch thrust-vector angle, tan 


f F ^ 

r N 


, deg 


nondimensional distance of the first grid 
normal to the surface 
turbulent energy dissipation 
laminar viscosity coefficient 
density, slug/ft^ 

Reynolds stress components 
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rOMPt ITATIONAL METHOD 
Governing Equations 

The PAB computer code solves the three- 
dimensional, Reynolds-averaged Navier-Stokes 
(RANS) equations and uses one of several turbulence 
models for closure of the RANS equations. The 
governing equations are written in generalized 
coordinates and in conservative form. In an effort to 
decrease computational resources, the simplified, thin- 
layer Navier-Stokes equations are implemented into 
PAB. This approximation neglects derivatives in the 
viscous terms streamwise and parallel to the surface, 
since they are typically negligible in comparison to 
the derivatives normal to the surface. Extensive details 
of PAB are found in references 28 and 29. 

The flow solver was written with three numerical 
schemes: the flux vector-splitting scheme of van Leer 
(ref. 30), the flux difference-splitting scheme of Roe 
(ref. 31), and a modified Roe scheme primarily used for 
space marching solutions. These schemes implement 
the finite volume principle to balance the fluxes across 
grid cells and the upwind biased scheme of van Leer or 
Roe to determine fluxes at the cell interfaces. Only 
the inviscid terms of the flux vectors are split and 
upwind differenced, while the diffusion terms of the 
Navier-Stokes equations are centrally differenced. The 
details and applications of these methods are given in 
references 29-3 1 . 

For this study and other typical three-dimensional 
simulations, the solutions are computed with the van 
Leer and Roe schemes. An iteration to steady state in 
a three-dimensional computational domain includes a 
forward and backward relaxation sweep in the 
streamwise direction, while implicitly updating each 
cross plane. 

Turbulence Modeling 

Turbulence modeling is required to predict 
solutions for many flow fields. The PAB code can 
perform several turbulence simulations by 
implementing either an algebraic or 2-equation, linear 
or nonlinear turbulence model. An algebraic 2-layer, 
Baldwin-Lomax model is accurate for simple viscous 
flows because the turbulent viscosity p T is determined 
by a local function. A 2-equation k-e model with 
second order closure is used to model more complex 
viscous flow features. A second equation is used to 
solve for the turbulent length scale in addition to the 
equation for turbulent kinetic energy (k). Since the k- 
e model has a singularity at solid surfaces, either a 
damping function or a wall function must be 
implemented to adjust the turbulent viscosity (£) near 


these surfaces. The grid in the boundary layer at wall 
surfaces must be well defined with a law-of- the- wall 
coordinate (y + ) of approximately 2 for adequate 
modeling of the boundary layer flow (ref. 32). The 
restriction on y + may be relaxed to 50 if a wall 
function is implemented. However, it is customary to 
restrict the use of wall functions to attached flows. 

Both linear and nonlinear turbulence simulations 
use the standard model coefficients of the k-e equations 
as a basis of formulation. The linear k-e turbulence 
model is an eddy viscosity model with the following 
formulation for v. 

T =T L + T r (1) 

a ij y ' 7 

where, 


T y — ^pL^kk^ij 2 Pl.S,j 

(2) 

rjj = ^(pk + p T S kk) S u ~ 2 Pt S „ 

(3) 

s 

9 2 ox, dx, 

(4) 

II 

( 5 ) 

The damping function, f M , is 

an empirical 


function, while is set to 0.09 for the standard linear 
k-E turbulence model. The turbulence model has one 
equation for turbulent kinetic energy, k, and one for 
turbulent energy dissipation, e. The k and £ pair ol 
coupled transport equations are written in conservative 
form which can be uncoupled from the Navier-Stokes 
equations and from each other to decrease 
computational requirements. In an effort to decrease 
numerical stiffness, the k and £ equations are solved at 
approximately 25 percent of the Courant-Friedrichs- 
Levy (CFL) number for the Navier-Stokes equations. 

The PAB code, with 2-equation, linear Reynolds 
stress modeling, was chosen for these simulations 
because this code was developed for and can accurately 
predict propulsive flows with mixing, separated flow 
regions, and jet shear layers. A modified Jones and 
Launder form (ref. 33) of the damping function (f M ) 
was utilized to treat the singularity at the wall. A 
high Reynolds number model with no damping 
function was implemented in the free stream blocks. 

Performance Calculation 

The PAB code contains a performance module (ref. 
34) that utilizes the momentum theorem applied to a 
user-defined control volume to calculate nozzle or 
aerodynamic performance. Quantities such as lift, 
drag, thrust, moments, heat transfer, and skin friction 
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may be computed for many complex geometric 
configurations and multi-stream flows. Each quantity 
is updated throughout the solution development to 
monitor convergence. 

Along flow-through sections of the control 
volume, mass and momentum fluxes, as well as 
pressure forces are integrated over each cell with 
equations 6 and 7. 

w=I(pU-N)AA (6) 

F flux = Z(pU (U-N) + (p - P~ )N}AA (7) 

where AA is the cell face area and N is the cell face 
unit vector. 

Along solid surfaces of the control volume, skin 
friction and pressure forces are determined. Surface 
pressure force F pressure is determined by multiplying cell 
static pressure by cell face area using equation 8. 

Fpressure = Z[(P-P~) N 1 M W 

The cell surface static pressure is calculated by 
extrapolating the cell centered static pressure to the 
surface where the velocity is assumed to be zero. 

The skin friction force F friction is calculated with 
only the velocity gradients normal to the surface 
contributing to the velocity terms of the viscous stress 
tensor. A two-point difference is used to determine a 
velocity gradient, one zero-magnitude velocity vector 
at the surface and a second at the cell center. 
Sutherland's formula (ref. 35) is used to calculate the 
dynamic viscosity at the surface by extrapolating the 
static temperature at a local cell center to the surface 
and using a reference viscosity and temperature 
condition. The total body force vector F is defined in 
equation 9. 

F —Fjj ux + F pressure + Ff riction (9) 

Ronndarv Conditions 

The PAB code has many options for defining the 
conditions of the inflow, outflow, free stream, wall, 
and centerline boundaries. For this study, Riemann 
invariants along the characteristics were implemented 
along the lateral and in-flow free stream boundaries. A 
first order extrapolation outflow condition was used at 
the downstream far field boundary. The primary nozzle 
flow and the secondary fluidic injection flow were 
specified with a fixed total temperature and total 
pressure boundary condition. A no-slip adiabatic wall 
boundary condition was implemented to obtain viscous 
solutions. 


Nozzle Geometry 

The two-dimensional, convergent-divergent 
(2DCD) fluidic thrust-vectoring nozzle was installed in 
an isolated nacelle to investigate the interaction 
between the external freestream flow and the nozzle 
exhaust flow. An isometric view of the isolated 
nacelle mounted on a strut is shown in figure 1 . The 
length of the isolated nacelle is 64.745 inches. The 
nozzle had a length of 8 inches, a throat area ot 
A =5.416 in 2 , an expansion ratio of A/A,=2.4, and an 
injection slot area of A #i#tor =0.565 in 2 . The design 
nozzle pressure ratio and exit Mach number, based on 
one-dimensional theory, are NPR D =14.588 and M=2.4, 
respectively. 

Computational Domain 

The computational mesh was fully three- 
dimensional with 2 blocks defining the internal nozzle, 

1 block representing the injection plenum and 16 
blocks representing the freestream domain. The 
internal blocking structure representing the nozzle and 
boattail is shown in Figure 2. The injection plenum, 
shown in Figure 3, had one-to-one grid matching with 
the nozzle divergent section mesh. A centerline cut ot 
the isolated nacelle and nozzle is shown in Figure 4. 
The far field was located 10 body lengths upstream and 
downstream of the aerodynamic nose and nozzle exit, 
respectively. The upper and lower lateral far field was 
located 6 body lengths above and below the 
aerodynamic body. The boundary layer was defined for 
a law-of-the-wall coordinate y + of 0.5 on the tine mesh 
spacing for adequate modeling of the boundary layer 
flow. 


RESULTS 

Three-dimensional computational solutions were 
predicted with the Navier-Stokes code PAB using 
turbulence closure and linear Reynolds stress modeling. 
Baseline, or no- vectoring simulations were computed at 
nozzle pressure ratio (NPR) and freestream Mach 
numbers listed in Table 1. Vectoring simulations were 
computed with static freestream conditions (Af=0.05) 
and with freestream Mach numbers of 0.3 to 1.2 at the 
NPRs listed in Table 1. In addition, two cases were 
computed at the design NPR D =14.588 with A/=0.3 and 
A/=1.2 freestream conditions. Vectoring simulations 
were computed with secondary to primary total pressure 
ratios (pjp t , p ) of 0.6 and 1 . The simulated conditions 
were chosen to coordinate with an experimental test 
matrix. 
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Baseline. No Ve ctoring Cases 

Four simulations with no vectoring (no secondary 
fluidic injection) were completed at the nozzle pressure 
ratios and corresponding scheduled freestream Mach 
numbers listed in Table 1. Internal performance is 
shown in Table 2. Discharge coefficient is dependent 
on the geometry in the convergent section and the 
geometry immediately downstream of the geometric 
minimum area (nozzle throat) for CD nozzles. 
Therefore, once the nozzle is choked, generally at 
NPR<2, discharge coefficient will reach a constant 
value. For the no vectoring cases, primary flow 
discharge coefficient ( C dp ) was relatively constant over 
the range of NPR. The boundary layer in the M= 0.7 
and NPR=5.2 case was slightly thicker than the other 
three cases, resulting in a 0.25 percent lower discharge 
coefficient. 

Thrust ratio (C f J increased nearly 8.5 percent as 
NPR increased from 3.6 to 7.2. At overexpanded 
conditions (NPR<NPR D ), losses due to shocks and 
separated flow inside the nozzle are significant. As the 
flow conditions approach the design NPR, the shock 
moves downstream in the nozzle and thrust approaches 
ideal. At fully expanded conditions, the nozzle is 
internally shock-free and losses from ideal thrust 
primarily include skin friction and flow angularity (due 
to nozzle divergence). 

There was a small pitch vector angle, <5 /? <0.4, for 
the cases with no secondary injection, see Table 2. The 
slot was included in the computational domain for the 
no injection simulations and left "open" to the primary 
nozzle. Therefore, small differences between the static 
pressures along the upper and lower wall existed (see 
figure 5(a)), and resulted in small non-zero pitch vector 
angles for the no injection cases. 

Mach contours and normalized internal static 
pressures along the nozzle centerline are shown in 
figures 5 through 8 for NPRs=3.6, 5.2, 6, and 7.2 at 
the scheduled Mach numbers, with no secondary 
injection. The pressure dropped in the convergent 
section of the nozzle (x<39.9 inches) as the subsonic 
flow accelerated to M= 1 at the throat (x=39.9 inches). 
The pressure increased at the throat as the flow 
negotiated through the sudden change in geometry from 
the convergent section to the divergent section of the 
nozzle (x>39.9 inches). As the flow accelerated to 
supersonic conditions in the divergent section of the 
nozzle, the pressure continued to drop. At these 
overexpanded conditions (NPR< 14.588), the pressure 
dropped below ambient pressure in the divergent 
section, but was adjusted to ambient conditions through 
a series of oblique shocks. At NPR=3.6 with M=0.3 


freestream flow conditions, the shock was located at 
x=42 inches, see figure 5(b). With no secondary fluidic 
injection, the pressures along the upper and lower wall 
were nearly identical, with the shocks located in the 
same axial location. The flow appears separated 
downstream of the shock, as evidenced by the relatively 
flat pressure distribution for x>42. 2 inches. The shock 
and flow separation is shown in the Mach contours in 
figure 5(c), with the shock located just upstream of the 
injection slot. The reversed flow at the point of 
separation from the upper wall is shown with u- velocity 
vectors in figure 5(d). As NPR was increased, the flow 
expanded further along the nozzle wall before the flow 
adjusted through a shock, compare figures 5 - 8. 

Hffcct of Freestream Flow on Pit ch Vectoring 

To understand the effects of the freestream flow 
on fluidic thrust vectoring, simulations were computed 
at nozzle pressure ratios of NPR=3.6, 5.2, 6, and 7.2 
with static freestream conditions (A7=0.05) and with the 
scheduled Mach numbers listed in Table 1. Internal 
performance for the simulations with static freestream 
conditions are shown in Table 3 and with the scheduled 
freestream Mach numbers are shown in Table 4. The 
ratio of secondary to primary total pressure was 0.6, 
resulting in a secondary injection flow rate of 5.8 
percent of the primary flow rate. 

The influence of the M=0.3 freestream flow acted 
to decrease vectoring effectiveness by increasing the 
negative pitch vector angle (S p ) 2.8 degrees at NPR=3.6, 
compared to static freestream conditions. Figure 9 
shows the internal, centerline pressure distributions for 
the NPR=3.6 cases with M= 0.05 and Af=0.3 freestream 
conditions. The shock on the upper wall moved further 
inside the nozzle, from x=42.5 inches to x=42.25 
inches, creating a larger region of separated flow with 
the influence of the freestream flow, figure 9(a). The 
freestream flow only influenced the pressure along the 
lower surface downstream of the injection slot, see 
figure 9(b). The pressure distributions along the wall 
downstream of the injection slot, were both at lower 
static pressures than ambient conditions. However, the 
jVf=0.3 freestream flow case had a greater suction force 
on the wall for x>42.2 inches, than the static freestream 
flow case. The combined effect of moving the shock 
upstream on the upper wall and increased suction on the 
lower wall resulted in a reduction of thrust vector 
effectiveness with a freestream flow of M= 0.3. 

Figures 10(a) and 10(b) show the Mach contours 
along the nozzle centerline for the NPR=3.6 cases at 
Af= 0.3 and M- 0.05, respectively. The larger region of 
separated flow along the upper wall with wind-on 
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conditions (A/=0.3 freestream flow) is evident in the 
Mach contours. In addition, the shock cells in the 
plume dissipated faster at M— 0.3 conditions than at 
static conditions. The Mach number reached A/=l .6 in 
the first shock cell of the external plume with static 
freestream conditions, but only reached A/ =1.4 with 
wind-on freestream conditions. 

The influence of M=0.7 freestream conditions also 
acted to decrease thrust vector effectiveness, with an 
increase in negative pitch vector angle of 2.9 degrees at 
NPR=5.2, compared to static freestream conditions. 
Figure 1 1 shows the effect of freestream flow on 
internal, centerline pressure distributions for the 
NPR=5.2 cases with A/=0.05 and M=0.1 freestream 
conditions. The freestream flow had the same effects on 
the pressure distributions that it had at NPR=3.6. The 
shock on the upper wall moved further inside the nozzle 
and a greater suction force was created on the lower 
wall, downstream of the injection slot with the 
influence of the freestream flow. The Mach contours 
for the NPR=5.2 cases are shown in figure 12. A shock 
formed and the flow separated from the upper wall with 
Af=0.7 freestream flow, figure 12(a), but the flow 
remained attached to the upper wall until the nozzle exit 
for the static freestream case. 

Trends similar to the NPR<6 cases are seen for 
the NPR=6 cases, see figures 13 and 14. However, at 
NPR=6, the suction on the lower wall, downstream of 
the injection slot is not as great as for the NPR<6 cases 
(figure 13). Therefore, the penalty in pitch thrust vector 
angle due to wind-on conditions is less than for the 
NPR<6 cases, 1.5 degrees compared to approximately 
2.9 degrees. 

The wind-on (Af=1.2) penalty to pitch thrust 
vector angle for the NPR=7.2 case was 2.2 degrees. 
Figure 15 shows the internal, centerline pressure 
distributions for the NPR=7.2 cases with M= 0.05 and 
A/= 1 .2 freestream conditions. There is little difference 
in pressures along the upper wall because the flow 
remained attached until the exit, figure 15(a). However, 
as with the NPR<6 cases, the suction on the lower wall 
downstream of the injection slot is much greater than 
the NPR=6 case, resulting in a greater penalty to thrust 
vector angle. The Mach contours along the nozzle 
centerline for the NPR=7.2 cases are shown in figure 
16. 

Figure 17 shows the effect of freestream flow on 
system thrust ratio, C ff(sys . Thrust ratio was decreased 
4.1 percent at NPR=3.6 and 0.83 percent at NPR=7.2 
with the influence of the freestream flow. As expected, 
the largest thrust penalty occurred at the most 
overexpanded condition. The further off design 
(overexpanded conditions) the nozzle, the more the 
internal flow separated from the divergent walls. The 


freestream flow influenced the separated, subsonic flow 
inside the nozzle, which pushed the shock further 
upstream and separated even more flow from the 
divergent walls. Although a thrust improvement is 
sometimes realized from shock induced separation, 
which effectively changes the expansion ratio to a more 
suitable one for a set of conditions, thrust was penalized 
in these cases. Loss of flow expansion along the walls, 
less flow momentum at the nozzle exit, and greater total 
pressure losses from shocks in the nozzle could explain 
the loss in thrust with the influence of the freestream 
flow. 

As expected, the freestream flow had minimal 
influence on internal nozzle discharge coefficient, see 
figure 18. Again, discharge coefficient is dependent on 
the geometry in the convergent section and the 
geometry immediately downstream of the throat. For 
these cases, the influence of the external flow remained 
in the divergent portion of the nozzle causing little 
impact on discharge coefficient. 

Effect of NPR on Pitch Vectoring 

Two simulations were completed to document the 
effect of NPR on pitch thrust vectoring at freestream 
Mach numbers of Af=0.3 and Af=1.2. The internal 
performance is shown in Table 5. Pitch thrust 
vectoring was achieved with a secondary to primary 
total pressure ratio of p,Jp tiP = 0.6. The secondary 
injection weight flow rate for p t Jp ttP = 0.6 was 5.8 
percent of the primary weight flow rate. 

At both Mach numbers, thrust ratio increased and 
pitch vector angle decreased with increasing NPR. In 
both cases, the shock on the upper wall moved outside 
of the nozzle from x=42.2 inches, while the shock and 
separation on the lower wall, upstream of the injection 
slot, also moved downstream as NPR increased, see 
figures 19 and 20. Eliminating the shock and 
separation along the upper wall improved thrust, but 
also decreased the pressure differential between the upper 
and lower surfaces that caused pitch thrust vectoring. 

Effect of Secondary to Primary Total Pressure Ratio on 
Pitch Vectoring 

The effect of secondary to primary total pressure 
ratio on internal performance is shown in Table 6. An 
additional simulation was computed at a nozzle pressure 
ratio of NPR=3.6 with Af=0.3 freestream conditions and 
a secondary to primary total pressure ratio of 1. The 
secondary injection weight flow rate was 9.65 percent of 
the primary weight flow rate. At NPR=3.6, the effect 
of increasing total pressure of the secondary injection 
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stream had an adverse impact on thrust vector angle and 
thrust ratio. 

In general, as injection flow rate increases from a 
no injection state, thrust vector angle also increases (ref. 
20, 23). However, at a particular injection rate, thrust 
vector angle will peak, resulting in diminishing returns 
for any further increase in injection rate. In some cases, 
adverse effects such as shock impingement (ref. 23) or 
flow reattachment will cause a decrease in vectoring 
with increased injection flow rate. 

In the case shown in Table 6 , pitch thrust 
vectoring effectiveness decreased as secondary to primary 
total pressure ratio increased from 0.6 to 1. The 
centerline pressure distributions shown in figure 21 help 
to explain the penalty to pitch thrust vector angle. The 
shock and flow separation from the upper wall moved 
further upstream, but then also reattached to the upper 
wall near x=42.8 inches, as injection total pressure 
increased. The flow reattachment to the upper surface is 
shown in the comparison of Mach contours in figure 
22. The shock in front of the injection slot on the 
lower wall moved further upstream as injection total 
pressure increased, see figure 21. The decreased pressure 
differential along the top and bottom wall resulted in a 
penalty of 3.38 degrees to pitch vector angle. 

System thrust ratio decreased about 1 percent with 
p t j Pt p = 0.6 (w,= 5.8 percent of w p ), but increased 

0.65 percent with p,Jp t ,p = 1 (w ? =9.65 percent ol w p ), 
compared to the no injection case (w v = 0.0 1 percent ol 
w ), see Table 6 . The flow reattachment to the upper 
wall in the pjp up = 1 case may have caused the thrust 
ratio improvement, but at a penalty to pitch vector 
angle. As expected, discharge coefficient was virtually 
unaffected by the increase of injection total pressure. 


CONCLUSIONS 

A computational study was completed on a fluidic 
pitch-vectoring, two-dimensional, convergent-divergent 
nozzle to document the effect of the external freestream 
flow on fluidic thrust vectoring effectiveness. The 
computational fluid dynamics code PAB was used with 
two-equation turbulence closure and linear Reynolds 
stress modeling. 

The external freestream flow decreased fluidic thrust 
vectoring effectiveness by 1.5 to 2.9 degrees over the 
range of NPR and Mach number schedule. The 
freestream flow influenced the separated, subsonic, 
internal flow by moving the location of the upper-wall 
shock upstream and caused increased suction on the 
lower wall, downstream of the injection slot. The 
overall decrease in pressure differential between the 


upper and lower walls resulted in less vectoring with 
wind-on conditions. 

Thrust ratio was decreased 4.1 percent at NPR=3.6 
and 0.83 percent at NPR=7.2 with the influence ot the 
freestream flow. As expected, the largest thrust penalty 
occurred at the most overexpanded condition where 
internal flow separation was the greatest. In addition, 
the freestream flow increased internal flow separation by 
moving the shock further upstream along the upper 
wall. As a result, thrust was reduced from a possible 
combination of reduced flow expansion, reduced flow 
momentum at the nozzle exit, and increased total 
pressure losses. 

The freestream flow had a minimal influence on 
nozzle discharge coefficient, as expected. The influence 
of the external flow was limited to subsonic flow in the 
divergent section of the nozzle, caused by separation at 
overexpanded conditions. Since discharge coefficient is 
dependent on the flow in the nozzle geometry up 
through the nozzle throat, discharge coefficient was not 
impacted by the external flow. 
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Mach NPR 8 P C fvsxx C (l C (ip C ds 

"5^5 3^6 - 1 9/16 .9048~ .9792 .9843 .898 1_ 

nnc cn -14.73 .9357 .9791 .9834 .9101 


.9785 .9827 I .9114 


.982 


Table 3. Internal performance for vectoring simulations 
with static freestream conditions, p t ,J P t,p = ®- 




Figure 2. Internal blocking structure 
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Table 4. Internal performance for vectoring simulations 
with the scheduled freestream Mach number, p,,JPt,p— 0.6. 

M NPR _Sp C_ d C d ,> 

-16.69 .8635 .9778 .9825 -9023_ 

~03~ 7.2 -9.98 .9413 .9770 .9810 .9132 


.588 


3 


.2 


14.588 


Table 5. Effect of NPR on internal performance for 
vectoring simulations, p, , =0.6. 


.9848 


Table 6. Effect of secondary fluidic injection on internal 
performance for NPR=3.6 with a freestream M= 0.3. 
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Figure 3. Injection plenum block. 
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Figure 4. Centerline <cut o^ the isolated nacelle and 
nozzle. 
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(a) Expanded view of the pressure distributions along 
the nozzle centerline, divergent walls. 



(b) Pressure distributions along the nozzle centerline. 



(c) Mach contours along the nozzle centerline. 



(d) Reverse flow on divergent wall, u-velocity vectors. 
Figure 5. NPR=3.6, A/=0.3, no injection. 




(a) 


Pressure distributions along the nozzle centerline. 


(b) Mach contours along the nozzle centerline. 
Figure 6. NPR=5.2, M= 0.7, no injection. 



(a) Pressure distributions along the nozzle centerline. 



(b) Mach contours along the nozzle centerline. 
Figure 7. NPR=6, M= 0.9, no injection. 
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(b) Mach contours along the nozzle centerline. 
Figure 8. NPR=7.2, Af=1.2, no injection. 




(b) Static (Af=0.05) freestream conditions. 
Figure 10. Mach contours along the nozzle centerline, 
NPR=3.6, p,Jp,„ =0-6 w/w p =5.8%. 
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(a) M=0.1 freestream conditions. 
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(b) Static (M= 0.05) freestream conditions. 
Figure 1 2. Mach contours along the nozzle centerline, 
NPR=5.2, pJp, tP =0.6 w/w p =5.8%. 



(a) M= 0.9 freestream conditions. 



(b) Static (Af=0.05) freestream conditions. 
Figure 14. Flow field comparisons along the nozzle 
centerline, NPR=6, pjp,, p =0.6 w/w„ =5.8%. 



(a) Upper wall 



(b) Lower wall 

Figure 13. Effect on freestream flow on internal, 
centerline pressure distributions, NPR=6, p,,/Pt, P —0.6. 



12 


American Institute of Aeronautics and Astronautics 


z, inches 










(b) Static (M=0.05) freestream conditions. 
Figure 16. Flow field comparisons along the nozzle 
centerline, NPR=7.2, pjp, p =0.6 w/w= 5.8%. 



Figure 17. Effect of a freestream flow on gross thrust 
ratio. 



Figure 18. Effect of a freestream flow on 
discharge coefficient. 
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(b) Lower wail 

Figure 19. Effect of NPR on internal, centerline 
pressure distributions with M= 0.3 freestream 
conditions. 
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Figure 20. Effect of NPR on internal, centerline 
pressure distributions with M=1 .2 freestream 
conditions. 



with A/=0.3 freestream conditions. 
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lecondary to primary total pressure ratio of 0.6. 

l_ . '1 



(b) Secondary to primary total pressure ratio of 1 . 

Figure 22. Effect of injection total pressure on Mach 
contours for NPR=3.6 with M=0.3 freestream conditions. 
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